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Abstract 

We consider a system of superconducting grains embedded in a normal metal. 
At zero temperature this system exhibits a quantum superconductor-normal 
metal phase transition. This transition can take place at arbitrary large con- 
ductance of the normal metal. 
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I. INTRODUCTION 



Quantum superconductor-insulator (or superconductor-normal metal) phase transitions 
can take place at zero temperature due to variation of parameters of the system. For example 
in experiments the transition takes place as a function of the degree of disorder in a 
superconducting film. In experiments of the transition was mediated by a magnetic 

field. 

It has been suggested that a disordered superconducting film can be described as a 
network of Josephson-coupled superconducting grains shunted by resistors P,r7|,PJTT[] . In 
this case the Coulomb interaction between electrons in grains suppresses fluctuations of the 
number of electrons in a grain and, due to uncertanty principle, increases the amplitude 
of fluctuations of the phase of the superconducting order parameter |]7|j9HTI|. The compe- 
tition between the charging energy and the Josephson inter-grain coupling energy leads to 
the phase transition. Shunting resistors play a double role in the model: a) dissipation in 
the resistors tends to suppress fluctuations of the phase of the order parameter b)tunneling 
between superconducting grains and the resistors renormalizes the capacitance of the grains 
and thereby the charging energy. As a result, in the framework of this model and in two 
dimensional case the onset of superconductivity corresponds to the normal state film's con- 
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ductance G of the film of order of ^. Since at G < ^ the system should be in the insu- 
lating state, it has been conjectured that the transition is of the superconductor-insulator 
type UsHTlj and that the behaviour of the system near the transition is universal. On the 
other hand, a renormalization group analysis |T^, which starts from a perturbation the- 
ory for slightly disordered uniform superconductors, showed in 2D case a zero temperature 
superconductor-metal transition. 



Recently Feigelman and Larkin |T3| reconsidered the problem in the framework of a 
model of superconducting grains embedded in a normal metal film. They showed that a) 
the superconductor-normal metal transition takes place and b) deep in the metallic phase 
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parameters of the system can be calculated with the help of perturbation theory. At G > ^, 
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however, the critical concentration of grains turned out to be exponentially small. 

On the other hand, the theories concerning granular superconductors PJ^,p|,p!I]] were 
done in the limit when the modulus of the order parameter on a grain does not fluctuate. In 
the case of no reflection on the superconductor-normal metal boundary this corresponds to 
R ^ ^. Here R is the grain's radius, ^ = mm[^; ^/^] is the zero temperature coherence 
length of the bulk superconductor, D = ^ is the electron diffusion coefficient in the metal, 
vp is the Fermi velocity, / is the electron elastic mean free path and Aq is the zero temperature 
value of the gap in the bulk superconductor. 

In this paper we consider the opposite case when R < ^ and show that the zero temper- 
ature quantum superconductor-metal transition as a function of grains' concentration can 
take place in samples with arbitrary large conductance and can exist even in the absence of 
any disorder in the sample and in the approximation when electrons in the normal metal do 
not interact. The critical concentration of the grains turns out to be relatively large. 

We will consider a quasi-two-dimensional film of thickness a ^ ^, which consists of 
superconducting grains of radius R embedded in a nonsuperconducting metal. The results 
can be easily generalized for d=3 case. We will assume that densities of states i' of the 
superconductor and the metal are the same and that the spatial dependence of the electron- 
electron interaction constant in the the sample has the form 

{Xs> if |r - rA < R 
(1) 
Aat < if |r - Til > R 

Here and Xn are electron interaction constants in the superconductor and in the metal 
respectively; r j are coordinates of centers of the superconducting grains and index i labels the 
grains. At zero temperature the linearized mean field equation Jl^ for the order parameter 
A(r) 

A(r-) = A(r-) J driK{r,ri)A{ri) (2) 

has a solution A(r) = f(\r — ri\) at R = R^-^\ Assuming no reflection at 
the superconductor-metal boundary we get R^-f^ ~ ^ [jl5|. Here K{r,ri) = 



/ deG{—e,r,ri)G{e,r,ri), G{r,ri) is the electron Green's function in the normal metal. 
It is convenient to normalize / (ir|/(r')p = 1. At i? < i?*,™-^-' the mean field value of the 
order parameter is zero (A(r) =0). To find the value of A(r) in the case R > R^^^'' one 



has to add to Eq.2 terms nonlinear in A ||T7|. 



II. QUANTUM FLUCTUATIONS OF THE ORDER PARAMETER IN AN 

INDIVIDUAL GRAIN 

We will consider the case R < R'^-^^ when for an individual grain Eq.2 does not have 
nonzero mean field solution. To calculate the correlation function of quantum fluctuations 
of the order parameter it is convenient to use a parametrization Aj(r,t) = ai{t)f{\r — rj|), 
which reflects the fact that, at {R'^^^ — R) <^ R'^^\ the amphtude of quantum fluctuations 
of ai(t) is large while the amplitude of fluctuations of the shape of /(r) is small. To describe 
dynamics of the order parameter in a grain we use the effective action 

'5. = ^ro/^H|a;| + -)|a.MP (3) 
J Ztt t 

whose derivation we outline in the Appendix. Here tq = min[:^] ^] is the the time of 
electron flight through the grain and 
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Rt^^ - R 



(4) 



Using Eq.3 we get 



< a.M,a:(-a;) >(o)= ^ i) (5) 



which in t representation corresponds to 



ift>r 

UTQ ^ t ' 



<a.(t),a*(0)>(o)=<| (6) 

Here the subscript "(0)" indicates that the correlation function is calculated for a sin- 
gle existing grain "i". Eqs.5,6 hold as long as terms nonlinear in |ajp in the effec- 
tive action can be neglected, i.e. if r ^ 5^^. Here 5 = {vRj^a)~^ is the average 



level spacing in the grain. They correspond to the casual Green's function Gc =< 
T(^^(r,t)^_^(r,t)^+(0, 0)^1^(0,0)) >, where a is the spin index. To get the retarded 
Green's function Gn =< 6'(t)[^+(r, t)^l^(r, t), ^^(0, 0)^_^(0, 0)]+ > one has to make an- 
alytical continuation of Eq.5 with respect to u and then to make a Fourier transform. As 
a result, Gji{t) ~ exp(— ;^). It is interesting that the asymptotic time dependence of Eq.6 
is the same as the one obtained in the case R ^ R^J^-^^ ||13| with the help of a complicated 
renormalization group analysis of the Caldeira-Leggett effective action In the latter 

case there is a non-zero mean field order parameter on a grain and the correlation func- 
tion decays with time due to phase fluctuations mediated by the interaction with quantum 
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electromagnetic fluctuations in conducting environment. At G > ^, however, the corre- 
sponding correlation time turns out to be exponentially large, which is different from the 
case R < R^""^^ Eq.4. 



III. THE SUPERCONDUCTOR-METAL TRANSITION IN A SYSTEM OF 

SUPERCONDUCTING GRAINS 

Let us consider now the case when the concentration of superconducting grains em- 
bedded into the quasi-two-dimensional metallic film is finite. To describe this system in the 
case when ctj are small we will use the effective action S = So + Smt where 5*0 = Y.iSi, while 
Sint describes inter-grain interaction via the normal metal and also the nonlinear in \ai\^ 
contributions to the action. 

Sint = - Yj J duJijai{u)a*{-uj) + j duidu2du'iai{uji)a*{uj2)ai{u^)a*{-uji - - tus) 

(7) 

where Jij have the meaning of Josephson coupling between grains i and j, and b ~ 

The mean field approximation corresponds to the minimum of 5* at a; = 0. If |rj — rj\ is 
small enough the solution of the Uzadel equation Eq.25 yields Josephson coupling between 
two grains of the form ( | — r j | ^ -R) 
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However, in the case of finite two-dimensional concentration such an expression for 
Jij would lead to a logarithmic divergence of the ground state energy density. 

On the other hand, at large |rj — rj\ electron- hole pairs diffusing through the metal 
between grains i and j will experience Andreev scattering from the superconducting grains 
situated between grains i and j. An example of such a grain "A;" is shown in Fig.l. Due to 
the Andreev nature of reflection |jl6| electrons scattered by the grain "/c" will be reflected 



into a hole moving in the direction opposite to the initial direction of electron's motion. 
Thus this electron will never reach the grain j. As a result we have 



"^J = W - vV ^"^^^ 



4"" = (-4^^)"' (10) 

Therefore the mean field equation for the ground state mean field value Q!(m/) of the order 
parameter on a grain has a form 

^ - 27rz//?2Arin(Ari/2L(-^)) + f,«2^^^ = q (11) 
Thus at T = an exponentially small mean field solution for the order parameter 

a(^;)~Ao/ai/2exp(--^) (12) 

exists at arbitrary small A^. 

However, at small enough A^ the amplitude of quantum fluctuations of the order param- 
eter Eqs.5,6 becomes larger than its mean field value. In this case the mean field theory is 
not applicable. To show that at small A^ quantum fluctuations destroy superconductivity 
completely and that the normal metal state is stable we use a perturbation theory procedure 
similar to the one used in |T3]. A requirement for the perturbation theory in therms of 



to be valid is the convergence of the integral 
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J < ai{t)a*{0) >dt< oo (13) 

which in our case follows from Eq.6. The integral in this case equals 

Again, the Josephson couplings of the form Eq.8 would lead to divergency of the pertur- 
bation theory. To cut off the divergence in the absence of the magnetic field one can consider 
the case |]13[ when there is a repulsion between electrons in the metal and Aat 7^ . Then 



in the two-dimensional case we have Jij ~ vE?r~'^[l + 2z/|Ajv|/n(-^)]~^ and the perturbation 
theory will converge on the lengthscale L\ ~ i?exp(^j^^). In the presence of a weak mag- 
netic field H the Josephson intergrain coupling decays exponentially on distances larger than 
the magnetic length Lh = \J^- Thus the cut off length relevant for the convergence of the 
perturbation theory with respect to the term in Eq.7 proportional to Jij is Lq = min[LA; Lh]- 
It gives small corrections to the Eqs.5, 6 as long as NR'^ln{LQN^^^) <^ ^. In the opposite 
case the ground state of the system is superconducting. Thus we can estimate a relation 
between the critical concentration of the grains Nc and their critical radius Rc from the 
equation 



I p(m/) p I 



(14) 



For example, in the case R < R^-^^ and at = we have the estimate 



Nc ~ ^^\Xn\ (15) 

This is different from the case R ^ i?^"*-^) where at G > ^ the critical concentration is 
exponentially small. The difference originates from the difference in correlation times. 

Neglecting the second term in Eq.7 and assuming for simplicity that superconducting 
grains form a square lattice we get an expression for the correlation functions 

< ak{uj)ai{-uj) >= —{~i\u\ + - + —Jij)k] 

= exp(~iqr) — — , -— — (16) 

where components of q are whole multiples of and L is the samle size. 



K{q) = E exp{iqr)K{ri) (17) 

At |q|iV"i/2 < 1 we have 

K{q) = -nDNln{{l + {qLof)LlN) (18) 

The mean field Eqs. 16-18 are vahd because the characteristic radius of interaction between 
grains N^^^'^ln{LQN^) ^ N^^^'^ is much larger than the average intergrain distance. 

Qualitatively the picture of quantum superconducting fiuctuations in the normal metal 
is similar to superconducting fluctuations in uniform metals at temperatures which are close 
to the critical one p|,|l3: 

a) Due to quantum fluctuations conductivity of the system is enhanced as compared to 
the normal metal value. It exhibits a big positive magnetoresistance. 

b) At small magnetic fleld the zero temperature Hall coefficient is suppressed as compared 
to its normal metal value |]22[. It also exhibits strong magnetic field dependence. 

c) Diamagnetic susceptibility is enhanced and exhibits a strong nonlinearity as a function 
of magnetic field. 

d) Energy dependence of the density of states at the Fermi surface has a dip, whose 
amplitude is magnetic field dependent. 

IV. CONCLUSION 

We have considered a model in which superconducting grains with the radius R < Rc are 
embedded in the normal metal and have shown that in this case there is a zero temperature 
quantum superconductor-normal metal phase transition as a function of and R. This 
transition is driven, primarily, by fiuctuations of the modulus of the superconducting order 
parameter. Though the parameters of the transition, in principle, depend on D, it exists 
even in the case when there is no disorder in the sample. The critical concentration of grains 
turns out to be relatively large. 
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Calculations presented above did not take into account localization effects pO|pT| . In 
d=3 case they are small. The question of whether in two-dimensional case the metallic state 
is localized requires additional investigation. On the other hand, as we have discussed, 
the transition can take place at very large sample conductances when the localization length 
is also very large. 

We would like to mention that in the presence of the electron-electron repulsion in the 
metal the metal-insulator transition can exist even in the mean field approximation described 
by the Eq.2. For example, in the case As ~ |AAr| and R < K^^ we have Nc ~ 

Finally, we would like to mention, that in our opinion experimental data on 2-D films 
do not contradict the possibility of a zero temperature superconductor-normal metal 
quantum phase transition. 



V. APPENDIX 

We start with a standard expression for the partition function in a superconductors (see 
for example P5|): 



D^D^*e-xY>{iS^ff) 



(19) 



S, 



eff 



drdt[ 



A(r) 



dr]Tr[A{r,t)Gr,irt,rt 



(20) 



where Gri{rt, r't') is a matrix Green function which is a solution of the Gorkov equation 

[-ij^ - e(^V)cr, + u{r) - r/A(r, f)]G',(r, t; r't') = (5(r - r')S{t - t')! (21) 



Here u{r) is the external potential. 



A 



A* 



(22) 



and / is a 2*2 unit matrix in Nambu space. We assume that u{r) has a white noise statistics 
with correlation functions << u{r) » and << u{r)u{r') >>= j^S{r — r'). Here brackets 
<<>> stand for averaging over realizations of the scattering potential. 
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Averaging Eq.21 over realizations of u{r), neglecting all weak localization and mesoscopic 
corrections, and making the diffusion approximation we get Uzadel equations for the normal 



Expanding Eqs. 20-24 with respect to A we get an expression for the effective action Eq.3 
for an individual grain. 

To get Eqs. 9, 10 we have to solve the Uzadel equation Eq.24 in the normal metal between 
the superconducting grains. This solution corresponds to an electron diffusion in the metal 
and Andreev reflections from the superconducting grains. In the case R ^ R^-^ we have 
|A(t)| = Ao and the Andreev scattering cross section equals to R^. Since in our case 
A{t) <^ Aq the cross section is of order of R^^^^ and we get Eq.9,10. 

This work was supported by Division of Material Sciences, U.S. National Science Founda- 
tion under Contract No. DMR-9205144. We would like to thank A. Andreev, S.Chakravarty, 
M.Feigelman, A.I.Larkin, D.Khmelnitskii, S.Kivelson, and F.Zhou for useful discussions. 



and the anomalous Green's functions (See, for example, [Q) 




(23) 



G2 + |F|2 = 1 



(24) 
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FIGURES 




FIG. 1. A system of superconducting grains embedded into a normal metal. S stands for the 
superconductor while N stands for the normal metal. Solid lines correspond to trajectories of 
electrons diffusing between grains i and j. The dashed line corresponds to the hole reflected from 
the grain k. 
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